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ABSTRACT 

An overview of some studies which have recently been carried out in [l]-[3] 
on the control of third -generation spacecraft, as modelled by the M^AT space 
vehicle configuration, is made. This spacecraft is highly non- symmetrical and 
lias appendages which cannot in general be assumed to be rigid. In particular, it 
is desired to design a controller for MSAT which stabilizes the system and satis- 
fies certain attitude control, shape control, and possibly station-keeping re- 
quirements; in addition, it is desired that the resultant controller should be 
robust and avoid any undesirable "spill-over effects". In addition, the control- 
ler obtained should have minimum complexity. 

The method of solution adopted to solve this class of problems is to formu- 
late the problem as a robust servomechanism problem [5] -[7], and thence to obtain 
existence conditions and a controller characterization to solve the problem. 

The final controller obtained for MSAT has a distributed control configura- 
tion and appears to be quite satisfactory. 

INTRODUCTION 

This paper summarizes studies carried out in [l]-[3] on control system 
structures known as third-generation spacecraft. Such spacecraft have: 

(1) Large mass 

(2) High power 

(5) Large non- symmetric flexible appendages 

(4) Precise communication RF bean control requirements. 

In particular, the class of spacecraft represented by the Mobile Communications 
Satellite (MSAT) is used as a reference for these studies. This spacecraft has 
non- symmetric appendages which cannot be assumed to be rigid (see Figure 1) . 

There are a number of control problems associated with the attitude-control, 
shape-control and possibly station-keeping control for such third generation 
spacecraft (referred to as LFSS) , which may be listed as follows: 

A . Ti.e LFS S Control Problem 

Problem 1: Lightly Lamped, Oscillatory Plan 1 ' 

A LFSS has eigenvalues either at the origin or approximately distributed a- 
long the imaginary axis. One of the basic objectives that a controller must ac- 
complish in this case is to stabilize the rigid body modes of the LFSS, and at 
the same time to stabilize the elastic modes of the LFSS. This is called the 
LFSS stabilization problem. 

•This work was supported by the Department of Communication? , Ottawa, Canada under 
contracts DOC-CR-SP-82-007, D0C-CR-SP-S3-002, DOC-CR-SP-84-002. 
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Problem 2: Modelling 

In modelling a LFSS, experience has shown that dynamic analysis may provide a 
framework for the model 1 ing of the low frequency elastic modes of the LFSS in a 
reasonably accurate way, but that the high frequency elastic modes cannot be ex- 
pected to be determined accurately, i.e. there will always be errors present in 
modelling the high frequency elastic modes of the LFSS. In addition, the calcul- 
ation of dampening effects on the LFSS can only be done with great uncertainty. 

Problem 3: The Infinite Dimensional Plant - The "Spill-Over Problem" 

The classical modelling of elastic structures as continua results in the well 
known "infinite dimensional" system reprc station of a LFSS. Whether or not one 
adopts this infinite dimensionality representation seriously from an engineering 
standpoint, there is no question that the number of system elastic modes present 
in a LFSS is always larger than the number which any design model of a LFSS can 
accommodate. In trying to control the modelled rigid and elastic modes, it is 
essential that the controller should not cause. these unmodelled high frequency 
elastic modes to become unstable. This is called the "Spill-Over Problem". 


Problem 4: The Sensor/Actuator Placement Problem 

The LFSS is intrinsically distributed, and the configuration of control hard- 
ware is not in general specified. Thus, unlike many conventional control problems, 
part of the LFSS control problem is in determining the number and location of 
sensor/actuators on the LFSS. 

Problem 5: Requirement for Multivariable Control Theory 

The concept of "third generation" spacecraft, unlike the first and some 
second generation spacecraft, precludes single-input, single-output control design. 
Some type of multivariable control design method is mandatory to deal with the 
severe interaction occurring in the system. 

Problem 6: Minimization of Number of Sensors/ Actuators 

Th. j is a practical limiation on the quantity of hardware that can be dis- 
tributed over the LFSS vehicle. This implies in particular that one cannot assume 
full state feedback is available, and that the number of actuators/ sensors used 
must be limited, i.e. one must minimize any unnecessary sensor/actuators required 
for LFSS control. 

The following problem definition is now given: 

B. The LFSS Robust Servomechanism Problem 

Assume that a LFSS can be exactly described by the following finite dimen- 
sional linear time invariant model: 

x = Ax + Bu + Eio 

y = Cx + Fto (1) 
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n <4 

where xeR is the state, ueR is the control (actuator inputs), y eR m are the 

r ® o 

measured (sensor) outputs, and yeR are the outputs to be regulated. Here u)eR 

are assumed to be constant unmeasurable disturbances applied to the structure, 
w * 

<^ m eR m are assumed to be constant unknown measurement errors and e = y-y ^ is the 

error in the system where y re £ is a constant set-point. Thus, it is assumed that 

(1) may include an arbitrarily large number of elastic modes (but not infinite). 

Assume now that an approximate model of (1), called the design model for (1), 
is given by: 


x = Ax + Bu + Eio 


y = Cx + Foj 


( 2 ) 


where xeR is the state of the design model, and where n«n. It is desired now 
to find a controller based on the design model (2), such that when it is applied 
to (1), the system is asymptotically stable, i.e. no spill-over occurs, and such 
that : 


lim e(t) = 0 , Vx(0)eR n , VtoeR , Vw eR m (3) 

m 

This is called the LFSS Robust Servomechanism Problem, which includes the follow- 
ing subproblems: 

(1) Stabilization 

(2) Station-keeping 

(3) Attitude control 

(4) Shape control. 


THE MSAT CONTROL PROBLEM 

The MSAT spacecraft is illustrated in Figure 1. It consists of four compon- 
ents, one of wich is rigid (the bus) and three of which are flexible (the solar 
array, the tower, and the reflector). The tower-reflector-hub hinge point is 
assumed to have a gimbal (see Figure 2) . 

The coordinates assumed for each of these substructures are as follows: 

(1) Bus - three rigid rotations (6 ,0 ,0 )• 

x y z 

(2) Tower - relative displacement of tower tip to tower root 

- relative angular displacement of reflector with respect to 
frame fixed at tower root (with zero gimbal angles) 

( a l ,a 2 ,C V ' 

(3) Reflector - two gimbal angles at tower-reflector-hub hinge point (8^8.,)' 

The actuators which are assumed to be available are as follows: 

(1) Eight thrusters f^, i=l,2,...,8, four from thrusters on the bus and four 

from thrusters at the reflector hinge point, aligned as shown in Figure 2. 
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( 2 ) 


Two torquers at the reflector hub, one about each gimbal axis 
Figure 2) . 




(see 


In this case, a design model and an evaluation model was developed in [4], 
in which the design model has 18 states consisting of S rigid body modes (corre- 
sponding to the three rigid rotations of the bus and two gimbal angles of the 
reflector) together with 4 elastic modes, and the evaluation model has 32 states 
consisting of 5 rigid body modes and 11 elastic modes. Table 1 gives the eigen- 
values of the open loop system for the two models. The models used in this study 
included the effect of dampening terms D, Dg (see Table 1) . 


TABLE 1 : Open Loop Eigenvalues of MSAT Vehicle 



Standard Design Model 

Evaluation Model 


With Damping 

With Damping 

With Damping 

With Damping 


Term D 

Term D 

Term Dg 

Term Dg 


Excluded 

Included 

Excluded 

Included 

Rigid 

Body 

Modes 

0 (repeated 

0 (repeated 

0 (repeated 

0 (repeated 

10 times) 

10 times) 

10 times) 

10 times) 


0±j0. 124 

-0.000923±j0. 124 

0±j0.124 

-0.000923±j 0.124 


0±j0. 239 

-0.00170 ±j0.240 

0±j0.151 

-0. 000 853±j 0.151 

Elastic 

0±j 0.556 

-0.00856 ±j0.556 

0±j0.239 

-0.00171 ±j0.239 

0+j 0.780 

-0.0211 ±j 0.779 

0±j 0.556 

-0.00856 ±j 0.556 

Body 



0±j0.690 

-0.00553 ±j0.690 

Modes 



0±j0.780 

-0.0211 ±j 0.780 



0±j 1.55 

-0.0751 ±jl.55 




0±j3. 14 

-0.0280 ±j3.14 




0±j3.96 

-0.0528 ±j3.96 




0±j9 .95 

-0.524 ±j 10. 1 




0±jl4.0 

-1.17 ±jl3.8 


It may be noted that the elastic modes of the evaluation model interweave with 
the elastic modes of the design model. 

A. Description of Problem to be Solved 

Ii, this case it is desired to solve the LFSS Robust Servomechanism Problem 
for the MSAT vehicle. In particular, there are two separate requirements for the 
controller to be designed for the MSAT vehicle: 

Requirement I 

Find a controller, based on the MSAT design model, which solves the following 
problems: 

• Stability: stabilize the 5 rigid body modes and the 4 elastic modes of the 

system. re f 

• Attitude control: regulate 0,9,0 to desired constant set points 0 , 

ref ref X y Z a 

0 , 0 respectively, in the presence of unknown constant disturb- 

y z 

ances. 
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• Shape control: regulate ®2 +a 2* ^ ^ 1 ’ ^ ^^2’ ^"^ 3 * a 3 to zer0 » zn 

the presence of unknown constant disturbances. 

• Spill-over problem: it is desired that the controller should satisfy the 

above requirements, and not cause any instability to occur with respect 
to any of the vehicle's elastic modes which are not included in the 
design model. 

• Controller complexity: it is desired to minimize the number of sensors and 

actuators which are required to solve the problem. 

• Discrete controller implementation: it is desired that the controller, when 

implemented digitally, should not require an excessively large sampling 
rate to maintain stability. 


Requirement II 

Apply the controller obtained, based on the MSAT design model, to the MSAT 
evaluation model, and verify that all objectives above are satisfied. 

The outputs to be regulated in this case are given by: 


y = ( 0 x *e ,e z ,6 1 +a 1 ,$ 2 +a 2 ,f" 1 6 1 ,f" 1 6 2 ,f" 1 6 3 ,a 3 )' 


(4) 


B. Assumptions Hade in Problem Formulation 

In this problem, it is assumed that there is no requirement for controlling 
the w x> uiy, rigid body modes. (Note: this assumption is not essential, e.g. 

12], 13] also deals with the case of station-keeping.) It is also assumed that 
there is no need to include any gyroscopic terms in the design and evaluation 
models . 


METHOD OF SOLUTION ADOPTED TO OBTAIN 
A CONTROLLER TO SOLVE PROBLEM 


The method of approach adopted to solve this problem was based on using the 
results of the "robust servomechanism problem" [5] -[7], in conjunction with a 
parameter optimization method [ 8 ] to determine the controller's parameters, e.v 
see [9] which solves a special case of the above problem when the sensors and 
actuators are collocated, using a decentralized control configuration. In this 
case, existence conditions for a solution to the problem were obtained, and a 
necessary controller structure developed. In particular, it was found that any 
controller which solves the MSAT problem specifications must consist of a "servo- 
cornensator" [5] (unique), together with a stabilizing compensator (non -unique ) . 
In this study, the simplest possible stabilizing compensator, i.e. a stabilizing 
compensator consisting of only proportional and rate feedback terms, was used. 

In this case, in order to satisfy the existence conditions obtained for a 
solution to exist to the problem, it was necessary to choose the following inputs 
(actuators) and measurable outputs (sensors) for the controller: 


Outputs (sensors) : 


y m = (e x ,6 y ,e z ,6 l’ e 2 ,0t r a 2 ,f 6 l* f V 


(S) 
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( 6 ) 


Inputs (actuators) : 

u * (8 c.' 8 e i - 8 6 2 - f i' f 2' f 5- f 6 ) ' 

* * * * * 

where g^ , fj, f 7 , f,., correspond to various combinations of the thrusters 
f j ,t' 2 , . .? ,fg (see Figure 2), as described in Appendix I. 


FINAL CONTROLLER CONFIGURATION OBTAINED 


In this case, »„e following distributed controller was obtained as a solu- 
tion to the MSAT robust servomechanism problem, based on the MSAT design model: 
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where s denotes the Laplace Transform operator, where 


fe 1 



f ref 
0 rer 

X 



X 

0 

4 

^J 2 

e ref 

y 

0 

Z J 


s+yj 

y 

e ref 

2 J 


and where Kj, K 2> Kj, K 4> K 5 , y are given as follows: 


1.43 

0.500 

24.7 

1.34 

-0.0460 

0.0255 

4.64 

1.12 

15.6 

-0.000439 

-6.81 

-0.000957 

0.00981 

-0.000483 

18.6 

0.00326 

38.0 

-0.231 

14.5 

-0.00955 

59.0 

-0.00916 

0.0216 

0.0127 

-2.40 
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K, 


K, = 


K 


4 


K 


5 


28.5 

0.510 

-136 

0.0653 

1180 

" 7. 14*10"' 
1.28*10*' 
-3.41*10*' 
1.63*10"* 
2.95*10"' 

-0.464 

-0.0144 

-0.00438 

0.136 

-0.000753 

0.0268 

0.0273 


10.0 

92.8 

-0.0191 

760 

-0.183 

2.50*10 

2.32*10 

-4.79*10 

1.90*10 

-4.58*10 

0.0433 

0.0536 

0.201 

0.0461 

1.13 

0.225 

-0.226 


0.0268 0.225 

0.0273 -0.226 


494 26.7 

22.3 312 

0.196 -0.00965 

-4.63 290 

0.432 0.254 


-0.920 

-0.00877 

372.2 

-0.191 

-48.1 


-4 

1 . 24*10~ 2 

6.68*10* 4 

-2.30*10* 5- 

-3 

5.58*10* 4 

7.80*10* 3 

-2 . 19*10* 7 

-7 

4.90*10" 6 

-2.41*10‘ 7 

9 . 31*10* 3 

-2 

-1.16*10* 4 

7.26*10" 3 

-4.78*10* 6 

-6 

1 .08*10" 5 

6.34*10~ 6 

-1 ,20*10” 3 _ 


Y = 2.0*10 


This controller is just a multivariable generalization of the classical 
three term controller used in classical control. The controller has minimal com- 
plexity ; . the sense that it has minimum order feedback dynamics and has the 
minimum number of actuators/sensors required in order to solve the problem. It 
is to be noted that no a priori assumption the distributed structure of (7) 
was made — the distributed structure of the controller (7) arose from the 
analysis automatically. 


PROPERTIES OF PROPOSED CONTROLLER 

The mam features of the proposed controller when applied to the MSAT design 
node! a evaluation model will now be described. The main features of interest 

are : 

(1) The stabilization properties of the proposed controller. 

(2) The steady state regulation px-operties of the proposed controller. 

The following results are obtained: 

A . Eigenvalues of Closed Loop System Using Proposed Controller 

Table 2 gives a listing of all eigenvalues obtained by applying the proposed 
controller (7) to the MSAT design model and evaluation models. 
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TABLE 2 : Listing of Closed Loop Eigenvalues Using Proposed Controller (7) When 

Applied to MSAT Design and Evaluation Models ' 


Standard Design Model 

Evaluation Model 

- 0 . 00047t j 0 . 0085 t 

-0.0024±j0.016 . * . . 

-0. 0051 ±j 0.022 body 

-0.0097±j 0.030 | es 

-0. P10± jO. 031 i 

-0. 0004 7tj 0.0085 t 

-0.0024± j0.016 . .* . 

-O.OOSlij 0.023 ri «^ body 

-0.0097±j 0.030 | es 

-O.OlOt j0.031 J 

-0.00014±j0. 124 1 

-0. 006 l±j 0.240 elastic 

-0.017tj0.557 

-0.029ij0.7S0 

: body 
les 

-0.00014±j0.124 

-0.00020±j0.151 

-0. 00611 j 0.240 

-0.017tj0.557 

-0. 00791 j 0.690 

-0. 00291 j 0.780 elaStl< 

-0. 1291 j 1.35 ■°‘ 

-0. 0671 j 3. 16 

-0. 0691 j 3. 95 

-S.51j8.88 

-O.Slljll .3 

body 

les 

-5.0*10"“ t 

-5.0*10"“ 

-5.0*10-“ 1 

c nxln -4 servo-compensator 

foxier* “°f s 

-5.0*10"“ 

-5.0*10"“ + 

-1.7*10" 3 t 

-5.0*10"“ 

-5.0*10"“ 1 

-5 0*10"“ servo-compensator 

-5.0*10““ “°f es 

-5.0*10'“ 

-5.0*10"“ ♦ 

-2.0*10“ 3 t 

-?'S* 1 lS" 3 feedforward 

-2’ 0*10" 3 controller 

-2 .0*10" 3 ® 0< ? es 

-2.0*10" 3 4- 

-2.0*10" 3 t 

-2.0*10~ 3 e * . 
2 0*10" 3 feedforward 

-2]o*10" 3 controller 

-2.0*10“ 3 “°f es 

-2.0*10" 3 4 


It is observed that the resultant closed loop system is asymptotically stable 

for both the design and evaluation models, i.e. no undesirable spill-over effects 

occur. It is also observed that the dominant time constant of the system is 

mainly associated with the servo-compensator modes. This implies that one would 

expect for the case of tracking, that the dominant time response of the system 

would be associated with the feedforward controller modes, i.e. TC, 4500 sec 

dora 

4 8 min., and for the case of disturbance rejection, that the dominant time of 

the system would be associated with the servo-compensator modes, i.e. TC . 4 

aom 

2000 sec 4 0.6 hrs. This result is verified in the simulation studies to follow. 
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B. Steady-State Values of Outputs Using Proposed Controller : Tracking Case 

Table 5 gives a summary of results obtained for the case of unit step func- 
tion tracking, when the proposed controller (7) is applied to the MSAT design and 
evaluation model. It is observed that all 9 outputs of the system are asympto- 
tically regulated to their correct values as desired. 

TABLE 3 : Steady-State Values of Outputs Using Proposed Controller (7) When 

Applied to Design and Evaluation Model - Tracking Case 



e ref =i 

X 

.y 

0 ref =l 

z 

0 

X 

1 

0 

0 

0 

y 

0 

1 

0 

0 

z 

0 

0 . 

1 

V“i 

0 

0 

0 

e 2 *a 

0 

0 

0 


0 

0 

0 


0 

0 

0 


0 

0 

0 

a 3 

0 

0 

0 


Note: Any | number | <10 is assumed to be zero. 


C. Steady-State Values of Outputs Using Proposed Controller: Disturbance 

Rejection Case 


Tables 4 and 5 give a summary of all results obtained for the case of dis- 
turbance rejection, when the proposed controller is applied to the MSAT design 
and evaluation models respectively. In this case, it is assumed that a unit step 
function change occurs for different disturbances corresponding to g c ,g c ,..., 


fg.fg defined in Table 6. It is observed that the first 7 outputs of the system 
are asymptotically regulated to zero, and that the remaining two outputs are 
approximately equal to zero in all cases, as is desired. . 


D. Sampling Rate Requirements for Digital Implementation of Proposed Controller 

If it is assumed that the proposed controller (7) is to be implemented 
digitally, then it is necessary that the sensor outputs and actuator signals be 
updated at a fast enough rate so as to guarantee closed loop stability, when the 
the controller is applied to the evaluation model. In this case, on assuming 
that the sensor and actuator signals are updated at the same rate, it was found 
that a sampling rate of at least 0.1 Hz must be used to implement the proposed 
controller. This requirement is not demanding. 
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Note: Any | number |< 10 ^ is assumed to be zero. 



Note: Any |number|<10 ^ is assumed to be zero. 
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SIMULATIONS OBTAINED USING PROPOSED CONTROLLER 
TO SOLVE MSAT PROBLEM 


This section gives some typical simulations of the closed loop system ob- 
tained by using the proposed controller (7) applied to the MSAT design and evalu- 
ation models. Additional simulation studies are given in [3]. 

TAf 

A. Example No. 1 (Attitude Control: 0^ =1) 

In this example, it is assumed that the system has zero initial conditions, 
that there are no disturbances present, and that a unit step function change of 

+1 occurs in the set point for 0 at t=0, i.e. 0 re ^=l, 0 re ^=O, 0 re ^=G. 

x x y z 

Figure 3 gives a plot of all 9 output variables y given by (4) when the con- 
troller is applied to both the design and evaluation model in this case. It is 
observed that the system’s response is almost decoupled, i.e. the output © x is 
approximately equal to its desired value of +l.at t i 50 min, and that all other S 
outputs are barely excited. 

Figure 4 gives a plot of the 7 control variables u given by (5) for this 
example. 

B. Example No. 2 (Disturbance Rejection: fj.=l) 

In this example, it is assumed that the system has zero initial conditions, 
that all set points are identically equal to zero, and that a unit step function 
change of +1 occurs at t=0 corresponding to a disturbance thrust fj=l, where f5 

is defined in Table 6. This example would correspond to a misaligned thruster 

associated with the proposed controller. 

Figure 5 gives a plot of all 9 output variables y when the controller is 
applied to both the design and evaluation model in this case. It is observed 
that the elastic modes of the vehicle are now excited, and that the output vari- 
ables are asymptotically regulated to zero in approximately 2.7 hours, which is 
consistent with the closed loop eigenvalues of the system given in Table 2. 

Figure 6 gives a plot of the 7 control variables u for this example. 

C. Example No. 3 (Disturbance Rejection: J g =l) 

This example is similar to Example No. 2 except that it is assumed that a 
unit step function of +1 occurs at t=0 corresponding to a disturbance thrust f g =l, 
where f„ is defined in Table 6. This disturbance is representative of an arbi- 
trary constant disturbance which may affect the system. 

Figure 7 gives a plot of all 9 output variables y when the controller is 
applied to both the design and evaluation models in this case. It is observed 
that the elastic modes of the vehicle are now also excited as they were in 
Example No. 2, and that the output variables are satisfactorily asymptotically 
regulated with the same time constant as in Example No. 2. 

Figure 8 gives a plot of the 7 control variables u for this example. 


11 



TABLE 6: Definition of Disturbances Assumed 



— 
Disturbance forces corresponding to the thrusters 
f r f 2’ f 5 ,f 6 ,f 0’ f 9 res P ectivel y of figure 2 


Disturbance torques corresponding to gg ,gg 
respectively about the gimbal axis 3^,3^ 2 

i c >g c >§ c 

n z 3 

Disturbance torques in the bus about the x,y,z 
axis respectively 


ROBUST PROPERTIES OF CONTROLLER DESIGN METHOD 

A study of the robustness properties of the proposed controller design method 
was carried out [3] . This was done by comparing the controller designs obtained 
using the proposed method to different design models of MSAT. It was concluded 
that the proposed design method appears to be quite insensitive to the type of 
design model used, e.g. all controllers obtained, when based on MSAT design models 
which had at least two dominant elastic body modes included, produced stable 
closed loop systems and give satisfactory tracking/regulation, when applied to 
the MSAT evaluation model. Other studies showed that the controller is robust 
with respect to evaluation models of arbitrary complexity. 

CONCLUSIONS 

This paper gives a brief summary of the work performed in [l]-[3]. In these 
studies, the control system design of a third-generation spacecraft, as modelled 
by the MSAT space configuration is studied. This spacecraft is highly non- 
symmetrical and has appendages which cannot, in gen3ral, be assumed to be rigid; 
the elasticity of these appendages makes the control system design particularly 
demanding. In particular, it is desired to design a controller for MSAT which 
stabilizes the system and satisfies certain attitude control , shape control ard 
possibly station-keeping requirements. In addition, it is desired that the resul- 
tant controller should be robust and avoid any "spill-over effects", i.e. it 
should satisfy the problems' specifications based on only an approximate design 
model for MSAT being available. In addition, the controller obtained should have 
minimum com pl exity , i.e. a minimum number of sensors/actuators should be used. 

The method of solution adopted to solve this class of problems was to formu- 
late the problem as a robust servomechanism problem and thence to obtain existence 
conditions and a controller characterization to solve the problem. In this case, 
the controller obtained must contain a servo-compensator together with a stabiliz- 
ing compensator. 

The final controller obtained for MSAT has a distributed control configura- 
tion, and appears to be quite satisfactory, i.e. extensive testing of the con- 
troller shows that the controller is indeed robust with respect to the choice of 
the design model, and that it satisfies all specifications of the problem state- 
ment. 
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APPENDIX I 


Definitions of g* ,fj,f*,f*,f* 

* 

g is defined in terms of thrusters f . , f_, f_, f. as follows: 
C 7 1 c 3 4 




* 



* 

f_ are definea in terms of thrusters f , , f„, f_, f. as follows: 
2 1 2 3 4 
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Figure 2: Assumed control inputs for MSAT spacecraft (taken from [4]), 
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